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EXECUTIVE  SUMMARY 


This  report  presents  calculations  to  explore  the  use  of  depolarization  of  scattered 
light  as  a  discriminator  of  biological  particles  in  an  aerosol.  The  discrete  dipole  approximation 
was  used  to  model  appropriate  shapes  of  several  different  types  of  particles. 

To  model  an  ensemble  of  spores,  we  used  experimental  parameters  appropriate  to 
Bacillus  cereus.  For  this  case,  two  different  but  related  assumptions  for  particle  shapes  are  shown 
to  give  different  and  distinguishable  graphs  for  the  depolarization  D(©,T)as  a  function  of 
wavelength  and  scattering  angle. 

The  use  of  two  different  length  distributions  to  simulate  scattering  from  two 
different  strains  of  Escherichia  coli  gave  rise  to  distinguishable  graphs  for  D(©,/l)  for  the  two 
strains.  Here,  the  same  shape  (hemispherically-capped  cylinder)  as  well  as  the  same  value  for 
other  parameters  (i.e.,  width  and  index  of  refraction)  was  assumed  for  the  two  strains. 

In  the  preceding  calculations,  we  averaged  over  random  orientations  and  size 
distributions  based  on  experimental  measurements  for  a  real  aerosol.  The  calculated  results 
predict  that  D(©,T)  would  provide  some  identification  of  aerosols  made  up  largely  of  single 
biological  particles. 

The  Mueller  matrix  elements  Sn(0)  and  S22(©)  are  identically  equal  for  a 
spherical  scatterer,  so  that  D(Q,X)  =  1  -  S22  /  Su  is  identically  zero  for  an  aerosol  composed  of 
perfect  spheres.  Spherical  background  particles  mixed  in  with  the  bacteria  affect  the  angular 
dependence  of  the  depolarization  D(@,T)only  by  affecting  the  average  Su(©)and  average 
522(0)  for  the  aerosol  in  the  same  way,  leading  to  a  predictably  small  effect  on  D(@,T)  in  this 
idealized  model  case.  The  scattering  predicted  for  such  a  mixture  probably  would  be  similar  to 
what  would  be  observed  from  spores  or  bacteria  in  an  aerosol  with  a  large  portion  of  liquid 
microdroplets. 


Clusters  of  particles  are  also  observed  in  biological  aerosols.  To  keep  the 
calculation  tractable  for  the  case  of  clusters,  individual  spores  making  up  each  cluster  were 
assumed  spherical,  and  calculations  of  D(0,T)  were  made  for  plausible  clusters  of  six  to  eight 
of  these  “spore  spheres.”  Results  are  presented  for  two  different  geometries  of  seven  spore- 
sphere  clusters  indicating  the  extent  to  which  the  tightness  of  the  cluster  affects  the  graph  of 
D(0,T). 
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PREFACE 


The  work  described  in  this  report  was  authorized  under  Contract  No.  W91ZLK- 
05-P-0691.  The  work  was  started  in  March  2005  and  completed  in  November  2007. 

The  use  of  either  trade  or  manufacturers’  names  in  this  report  does  not  constitute 
an  official  endorsement  of  any  commercial  products.  This  report  may  not  be  cited  for  purposes 
of  advertisement. 

This  report  has  been  approved  for  public  release.  Registered  users  should  request 
additional  copies  from  the  Defense  Technical  Information  Center;  unregistered  users  should 
direct  such  requests  to  the  National  Technical  Information  Service. 
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CALCULATIONS  OF  LIGHT  SCATTERING  MEASUREMENTS 
PREDICTING  SENSITIVITY  OF  DEPOLARIZATION 
TO  SHAPE  CHANGES  OF  SPORES  AND  BACTERIA 


1.  INTRODUCTION 

There  has  been  practical  interest  in  determining  whether  polarized  light  scattering 
would  provide  important  input  into  determining  the  nature  of  an  aerosol  containing  biological 
particles.  This  is  of  interest  both  for  lidar  backscatter  and  for  point  measurements  where  samples 
of  the  aerosol  are  drawn  through  an  instrument. 

Bacterial  scatterers  of  light  in  an  aerosol  differ  in  size,  refractive  index,  or  shape 
from  the  background  particles  that  may  be  present.  Many  of  the  background  scatterers,  for 
example,  are  flat.  Other  scatterers  are  irregular  in  shape  but  are  often  modeled  in  calculations  as 
spheres.  Various  particles  may  also  clump  together  to  form  clusters.  However,  bacterial  particles 
are  ellipsoidal,  or  have  the  shape  of  a  prolate  spheroid,  or  often  have  a  shape  well  approximated 
by  capped  cylinders.  Depolarization  in  scattering  may  help  classify  aerosol  particles  according  to 
their  shape.  For  example,  scattering  from  perfect  spheres  produces  no  depolarization,  while 
scattering  from  ellipsoids  and  many  other  non-spherical  shapes  can  produce  considerable 
depolarization.  Thus,  the  polarization  of  the  backscattering  has  particular  interest  as  a  potential 
probe. 


Recent  experimental  work  examined  the  polarization  of  backscatter  signals  for 
aerosolized  BG  ( Bacillus  otrophaeus)  spores  mixed  with  naturally  occurring  organic  and 
inorganic  background  aerosol  materials.1  Certain  polarization  ratios  of  the  backscattered  light 
appeared  to  be  potentially  useful  as  discriminants.  A  particularly  interesting  variation  in  the 
energy  normalized  depolarization  ratio  at  1551  nm  divided  by  the  depolarization  ratio  at 
3389  nm  was  observed  and  might  be  useful  in  identifying  aerosols  containing  bacterial 
endospores.  Moreover,  it  was  proposed  that  the  large  imaginary  part  of  the  refractive  index 
around  3389  nm  might  suppress  depolarization  by  the  bacterial  component  near  that 
wavelength.1 


In  the  present  work,  we  examine  similar  questions  using  a  theoretical- 
computational  approach  based  on  modeling  the  scattering  of  aerosols  containing 
microorganisms,  including  mixtures  of  microorganisms  and  background  particles.  Such  an 
approach  based  on  computational  modeling  offers  the  advantage  of  allowing  us  to  study  the 
variation  of  the  scattering  with  parameters  in  ways  that  may  not  be  feasible  experimentally  (e.g., 
by  computationally  or  analytically  varying  the  ratios  of  the  different  particulates  once  the 
scattering  patterns  for  each  has  been  determined). 

Our  basic  approach  is  to  compute  the  set  of  Mueller  matrix  elements  {s ' }  for 

scattering  angles  between  0°  and  1 80°,  for  each  model  particle  type  /  of  interest,  averaging  over 
orientations  and  where  appropriate  over  size  distributions.  Then,  with  /?/  denoting  the  fraction  of 
particles  of  type  /,  the  appropriately  averaged  Mueller  matrix  element,  Sy,  can  be  obtained  as 
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(1) 


(s„)=Im; 

/ 

where  . S  j  for  each  species  includes  an  average  over  particle  orientations. 

The  methods  we  use  for  the  non-spherical  particles  are  an  extension  of  those  in 
previous  studies  by  the  present  authors  and  co-workers  who  applied  the  discrete  dipole 
approximation  (DDA),  also  known  as  the  "coupled-dipole  model,"  to  calculate  scattering  matrix 
elements  for  light  scattering  from  randomly  oriented  targets.2’3’4  For  the  simple  case  of  spherical 
background  scatterers  in  some  of  our  model  calculations,  we  use  the  Mie  scattering  program  of 
Bohren  and  Huffman.5 

Our  previous  studies  used  the  DDA  to  examine  light  scattering  patterns  from 
randomly  oriented  Escherichia  coli  bacteria  in  vivo  by  modeling  them  as  capped  cylinders  with  a 
distribution  of  lengths  and  diameters.  Detailed  comparison  with  experimental  data  gave  good 
agreement  and  supported  the  conclusion  that  observation  of  the  ratio  of  light  scattering  (Mueller) 
matrix  elements  <S34>/<Sn>  provides  a  rapid  and  convenient  method  of  monitoring  bacterial 
diameters  in  vivo  even  when  there  is  an  appreciable  distribution  of  lengths  in  the  bacterial 
population.  In  many  respects,  the  present  modeling  is  an  extension  of  this  previous  work;  in 
particular,  our  model  calculations  of  the  depolarization  from  E.  coli  based  on  experimental  size 
distributions  in  Section  4.2  builds  upon  the  earlier  work. 

The  biological  material  we  considered  in  the  present  work  included  spheroidal 
spore  shapes,  bacteria  (e.g.,  E.  coli,  whose  shape  can  be  approximated  as  that  of  a  capped 
cylinder),  and  spore  clusters.  Because  realistic  background  particles  occur  in  various  sizes,  we 
performed  model  simulations  with  experimentally  observed  size  distributions  and  random 
orientations  of  bacteria  and  idealized  background  particles. 

In  Section  2,  we  present  some  of  the  relevant  theoretical  background. 


2.  THEORETICAL  CONSIDERATIONS 

We  assume  light  incident  along  the  z  axis  with  the  coordinate  system  defined  in 
accordance  with  the  conventions  given  in  Bohren  and  Huffman5  and  illustrated  in  Figure  1 .  The 
scattering  plane  is  the  plane  containing  the  z  axis  and  the  scattering  direction  vector  labeled  k. 

The  polarization  state  of  either  the  incident  or  scattered  light  is  described  by  the 
Stokes  parameters  /,  Q ,  U ,  and  V.  The  first  parameter  /  is  the  total  irradiance  that  would  be 
observed  without  a  polarizer.  The  second  parameter  Q  is  the  difference  in  irradiances  polarized 
parallel  to  the  scattering  plane  and  perpendicular  to  the  scattering  plane,  and  would  be  measured 
using  parallel  and  perpendicular  polarizers.  The  third  parameter  U  is  obtained  by  rotating  a 
parallel  polarizer  by  +45°  and  -45°  and  calculating  the  difference  in  measured  irradiances.  The 
fourth  parameter  V  is  the  difference  in  right  and  left  circularly  polarized  irradiances. 
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Figure  1.  Scattering  coordinate  system.  The  light  is  incident  along  the  z  direction,  and 
the  scattering  is  calculated  in  the  direction  k,  with  scattering  angle  0 


2.1  Scattering  Matrix  Elements 

A  calculation  of  the  scattered  fields  for  scattering  in  a  particular  direction  of 
incident  light,  in  particular  polarization  states,  leads  to  a  linear  relation  between  the  Stokes 
parameters  of  the  incident  and  scattered  light.  The  specification  of  the  coefficients  in  this  relation 
describes  the  scattering.  For  example,  the  relation  between  field  polarizations  amplitudes  En 

parallel  to  the  scattering  plane  and  E±  perpendicular  to  the  scattering  plane  for  input  and 
scattered  light  are  given  by  the  amplitude  scattering  matrix 

The  electric  field  vector  on  the  right  hand  side  of  eq.  2  refers  to  the  incoming  light,  and  that  on 
the  left  is  the  field  for  the  scattered  light. 

Each  of  the  elements  of  the  amplitude  matrix  is  complex  valued.  The  phase  of  the 
matrix  elements  almost  never  can  be  determined  relative  to  other  particles  in  an  ensemble. 
Therefore,  we  confine  our  attention  to  calculations  that  do  not  depend  on  the  relative  phases 
(e.g.,  by  disregarding  multiple  scattering).  Single  particle  scattering  dominates  measurements  in 
all  observations  of  aerosols  where  the  particle  density  is  small,  as  is  the  usual  case. 


4  'l  ( S, 

v4i ) 


-  ikr 
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The  4x4  Mueller  matrix  for  a  single  scatterer  contains  all  the  polarization  and 
amplitude  information  of  the  light  scattered  from  either  a  single  particle  or  an  ensemble  of 
particles  that  receive  the  same  input  light.  The  elements  of  that  matrix  may  be  derived  from  the 
elements  of  the  amplitude  matrix  for  a  single  particle.  The  simple  equations  relating  the  elements 
of  these  matrices  are  given  in  Reference  5.  All  elements  of  the  Mueller  matrix  are  real,  and  are 
additive  for  all  the  scattering  particles.  Thus,  we  may  calculate  the  Mueller  matrix  for  each 
particle  in  an  ensemble  and  then  average  over  all  of  these  matrices  to  get  the  Mueller  matrix  for 
the  ensemble. 


The  Mueller  matrix  elements  {Sy}  describing  the  scattering  in  a  particular 
scattering  direction  (0,  are  then  the  1 6  elements  of  the  matrix  relating  the  Stokes  parameters 
for  the  incoming  light  and  the  scattered  light  in  the  scattering  direction 
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If  the  values  of  the  fields  in  eq.  2  are  evaluated  computationally  for  each  scattering  angle  of 
interest,  then  eq.  2  allows  the  four  scattering  amplitudes  Sj  to  be  calculated  at  each  scattering 
angle,  and  the  Mueller  matrix  elements  follow  from  the  four  Sj  amplitudes.  For  example,  Sn 
measures  the  total  scattering  for  unpolarized  incident  light 

Sn-itar+N’+N’+N2).  <4> 


while  the  value  of  S22  depends  on  changes  in  polarization 

sa=i(|s,f +N!-I5jf-N2)-  <s) 

In  particular,  the  depolarization  ratio  D,  defined  as 

D  =  (6) 

is  of  special  interest  in  the  present  work.  Note  that  S22  differs  from  Su  only  because  of  the 
presence  of  nonzero  values  of  |S3|"  +  |S4|’  according  to  eqs.  3-5.  But,  S3  and  S4  are  seen  in 

eq.  2  to  describe  the  change  from  polarization  in  one  direction  to  the  polarization  in  the  cross 
direction  to  it.  For  perfectly  spherical  scatterers,  S3  =  S4  =  0,  which  gives  no  change  in 
polarization  direction.  This  means  that  S22  =  Sn  for  such  scatterers.  Hence,  the  depolarization 
ratio  D(0,/l)  is  identically  equal  to  zero  at  all  scattering  angles  0  when  the  scatterers  have  an 
exact  spherical  symmetry. 
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2.2 


Coupled-Dipole  Method 


The  Mueller  matrix  elements  are  computed  by  calculating  the  scattered  fields  for 
the  various  polarization  of  incident  field  in  eq.  3.  This  determines  the  scattering  amplitudes,  and 
therefore  the  Mueller  matrix  elements  of  interest.  For  the  case  of  spherical  scatterers,  the  present 
work  uses  the  Mie  solution  to  calculate  the  scattered  fields  and  Mueller  matrix  elements  for  each 
scattering  angle.  For  other  particle  shapes,  we  use  variants  of  the  coupled-dipole  method  first 
introduced  by  Purcell  and  Pennypacker,2  and  also  known  as  the  DDA  method. 

The  DDA  method  is  based  on  replacing  the  scatterer  by  a  finite  array  of 
polarizable  sites,  which  can  be  conveniently  placed  on  a  cubic  lattice  within  the  volume  of  the 
scatterer.  The  polarizable  sites  are  subject  to  an  incident  oscillating  electromagnetic  field.  Let  E0 
be  the  vector  amplitude  of  the  incident  electric  field  E(r,/)  =  E0  exp[/(k  •  r  -  cot)]  where  k  is  the 

wave  vector,  and  the  location  of  dipole  /  is  at  position  r,.  The  electric  field  that  the  dipole  p ;  at 
position  Tj  produces  at  site  /  located  at  position  r,  is 

E(r,)  =  Vp,  (7) 

where 

p,  =«E(r.)  (8) 

and  where  a  is  the  polarizability  for  the  dipoles.  The  resulting  field  at  r*  from  all  the  other  dipoles 
and  from  the  applied  (complex)  field  E0(r, )  is 

E(r, )  =  Eo  (r, )  +  X  T</  *  aE(rj  )  •  (9) 

j 

In  terms  of  the  components  involved. 


m)  =  £„<'£)  +  ±±aT(i4;jv)-EUv) . 

7=1  r=l 


(10) 


where  T  denotes  the  complex  dipole-dipole  coupling  matrix  for  the  N  dipoles,  for  wavelength  A, 
and  refractive  index  n  of  the  particles.  Note  that  T  depends  on  wavelength,  and  on  the  relative 
locations  of  the  dipole  sites,  but  does  not  depend  on  the  refractive  index  n  or  on  the  polarizability 
a  derived  from  n  via  the  Clausius-Mossotti  relation.  Specifically,2 


m;jv)  = 


exp  (ikr9 ) 
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k'ri/ 


ro 
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3 (ikr0  -l) 
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(ii) 


A  matrix  element  of  T  is  labeled  here  as  T(/£yr)  corresponding  to  interaction  between  the 
dipole  i  and  the  dipole  j  (for  £=  1,2,3  labeling  ,  respectively,  the  x,  and  r  coordinates,  and 
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similarly  for  v  =  1,2,3 ).  Also  (r^  denotes  the  £  component  of  the  displacement  vector  of  site  / 
relative  to  site  j.  Equation  (10)  can  be  rewritten  as 

Z Z^A  -aT(i^Ju)\  E{jo)  =  £<,(/£)  (12) 

7=1  «=1 

which  is  the  equation  that  is  solved  to  find  the  field  and  the  dipole  moment  at  each  site.  It  can  be 
expressed  in  matrix  form  as 

|1-«T|.E  =  E0  (13) 


with  1  denoting  the  unit  matrix. 

We  assume  a  to  be  determined  by  the  refractive  index  n  of  the  continuous 
dielectric  particle  being  represented  through  the  Clausius-Mossotti  relation 


3  a2  - 1 
Mn  n~  +  2 


(14) 


where  M  is  the  number  of  dipole  sites  per  unit  volume. 

Although  various  refinements  and  improvements  of  the  DDA  method  have  been 
proposed,  the  basic  DDA  procedure  just  described  has  led  to  reasonably  accurate  calculation  of 
Mueller  matrix  elements  using  surprisingly  small  numbers  of  dipoles  to  represent  the  scatterer, 
and  has  proven  useful  for  modeling  scattering  from  bacterial  systems.3 


3.  PARAMETERS  AND  ASSUMPTIONS 

In  this  section,  we  describe  some  model  biological  particles  used  in  this  report, 
and  present  the  rationale  for  the  optical  parameters  and  distributions  we  use.  The  models  and  the 
parameters  were  chosen  as  examples  to  explore  and  illustrate  the  possible  importance  of 
depolarization  measurements  as  a  means  to  optically  characterize  some  representative 
bioparticles  and  backgrounds,  and  are  not  meant  to  be  exhaustive.  In  this  project,  we  examined 
the  behavior  of  certain  particular  Mueller  matrix  elements,  and  especially  the  depolarization 
ratio,  D(0,/l). 


In  many  biological  aerosols,  single  particles  dominate.  Therefore,  we  looked  at 
light  scattering  from  ensembles  of  spores  modeled  as  single  particles,  and  from  some  examples 
of  aerosols  of  single  vegetative  bacteria.  We  consider  Bacillus  endospores  (hereafter  simply 
called  spores).  There  are  other  types  of  spores  and  other  types  of  bacterial  spores,  but  these  are 
not  as  important  right  now  for  our  purposes.  The  Anthrax  and  other  Bacillus  spores  are  more 
precisely  named  endospores,  but  are  often  just  called  spores. 
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3.1 


Bacillus  Endospores 


For  the  spores,  we  used  parameters  appropriate  to  Bacillus  cereus.  The  size 
distribution  for  these  spores  is  excerpted  from  a  large  unpublished  set  of  electron  microscope 
measurements  due  to  Z.  Z.  Li,  J.  Czege,  and  B.  V.  Bronk.  These  measurements  are  shown  in 
Figure  2. 
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Figure  2.  Length  (left)  and  diameter  distribution  for  Bacillus  cereus  spores. 

The  smoothed  distributions  under  the  heavy  fitted  log-normal  graph  are 
representative  for  these  spores.  Inspection  of  the  microscope  pictures  indicated  that  the  small 
material  not  associated  with  the  fitted  curve  was  from  broken  spores  or  other  debris.  The 
diameters  (which  are  the  widths  in  Figure  1)  and  the  lengths  were  approximately  correlated  (i.e., 
longer  spores  generally  have  larger  diameter).  In  obtaining  the  graphs  of  Figure  2,  a  small 
positive  correlation  was  found  between  lengths  and  diameters  of  the  spores.  To  simplify  our 
calculations  and  reduce  computation  time,  we  assumed  the  correlation  was  exact  (i.e.,  with  a 
correlation  coefficient  of  1 .0).  With  this  simplifying  assumption,  we  derived  five  bar  histograms 
from  the  fitted  graphs  to  be  used  when  size  averaging  the  Mueller  matrix  elements  over  single 
spores.  These  histograms  are  shown  in  Figure  3.  In  our  averaging  over  sizes,  the  leftmost 
diameter  is  used  with  the  leftmost  length  and  so  on. 

In  addition,  we  needed  real  and  imaginary  parts  of  the  spore  index  of  refraction  to 
allow  us  to  do  calculations  that  should  provide  results  typical  for  B.  cereus  or  other  anthrax-like 
spores.  For  this,  we  used  interpolated  values  from  the  unpublished  experimental  data  of 
M.  Querry  (written  communication  from  M.  Milham).  These  values  are  given  in  Table  1. 

In  this  study,  we  use  the  above  parameters  for  B.  cereus  spores  to  see  what 
differences  in  D(@,i)  occur  for  two  different  plausible  shapes  of  an  individual  spore. 

In  electron  microscope  pictures  of  these  spores,  the  individual  spores  appear  like 
cylinders  with  rounded  ends.  Therefore,  in  our  first  set  of  calculations,  we  used  right  circular 
cylinders  with  radius  a.  The  model  cylinders  have  hemispherical  end-caps  of  the  same  radius  a 
where  2  a  corresponds  to  the  diameters  in  Figure  3.  If  the  length  of  the  cylindrical  part  is  (,  the 
overall  length  is  then  L  =  t+2a  where  this  quantity  assumes  the  value  indicated  in  Figure  3  for 
the  lengths. 
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Figure  3.  Distribution  of  diameters  (left,  black)  and  lengths  used  in  modeling  B.  cereus  spores. 
The  simplifying  assumption  of  direct  correlation  of  length  and  width  (smallest  length 
with  smallest  width,  etc.)  is  used  in  the  model. 


Table  1.  Refractive  Index  of  B.  cereus  Spores 
at  Representative  Wavelengths 


Wavelength  (nm) 

Real  (n) 

Im  (n) 

266 

1.619 

0.0599 

355 

1.570 

0.000969 

410 

1.54872 

0.000524 

1551 

1.503 

0.0004 

3389 

1.511 

0.0133 

In  a  second  set  of  calculations,  we  modeled  the  spores  using  the  shape  of  a  prolate 
spheroid.  For  the  prolate  spheroid,  the  defining  equation  is  given  by 

(x/o)2  +{yla\  +(z/b)2  =1;  wherea</>  (15) 

and  b  is  the  semi-major  axis.  Then  2b  is  the  length  corresponding  to  L  for  the  capped  cylinders, 
while  a  is  the  semi-minor  axis  and  the  radius  of  the  hemispherical  end-cap;  hence,  2  a  was  set 
equal  to  the  values  of  the  diameters  in  Figure  3. 
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In  each  of  these  two  cases,  we  used  the  DDA,  with  single  dipole  locations  filling 
the  space  inside  the  outer  boundary  defined  by  the  appropriate  equation.  Roughly  3000  dipoles 
per  spore  were  used  with  their  location  generated  as  described  in  Reference  2. 

3.2  Vegetative  Bacteria 

Vegetative  bacteria  come  in  many  forms.  However,  herein,  we  consider  just  one 
of  the  important  shapes.  That  is  rod-shaped  with  end-caps  similar  to  hemispheres.  In  earlier 
work,  we  used  microscopically  measured  diameters  and  lengths  as  inputs  for  Mueller  matrix 
calculations  based  on  the  DDA  model  and  the  capped-cylinder  shape.  The  graph  we  obtained  for 
S34  /  Su  gave  a  good  fit  to  the  corresponding  graph  obtained  experimentally.3 

For  the  present  calculations,  we  used  parameters  suggested  by  those  measured  for 
two  different  strains  of  Escherichia  coli.  These  and  other  vegetative  bacteria  have  decidedly 
different  refractive  indices  from  Bacillus  spores.  A  refractive  index  n  =  1.373  +  0.000396/  was 
measured  previously3  and  is  a  good  approximate  value  through  the  range  we  model  (355  nm  or 
400  nm  wavelength  of  scattering  light)  although  a  slightly  larger  imaginary  part  (0.00097)  was 
used  in  these  calculations  for  vegetative  bacteria. 

In  our  first  set  of  calculations  relating  to  these  cells,  we  averaged  over  orientation 
as  well  as  the  previously  measured  length  distributions  for  log-phase  cells  of  the  bacterial  strains 
B/r  and  K12  of  Escherichia  coli,  but  arbitrarily  used  only  a  single  diameter.  The  five  bar 
histograms  shown  in  Figures  4  and  5  represent  the  log-phase  lengths  measured  previously3-4  for 
the  two  different  strains.  We  used  these  log  phase  distributions  as  a  model  for  what  might  be 
encountered  in  an  aerosol,  and  we  used  the  distribution  to  average  over  sizes  in  our  calculations. 
The  diameters  also  have  a  distribution;  however,  they  were  uncorrelated  with  length,  probably 
because  these  rod-shaped  cells  grow  by  elongation.  The  K12  bacteria  were  found  to  be  wider 
than  the  B/r  bacteria  (averaging  about  1 .05  pm  for  K12  vs.  -0.85  pm  for  B/r).  However,  in  the 
present  calculations,  we  were  trying  to  see  how  changing  just  length  would  affect  the 
depolarization.  Therefore,  we  used  the  same  diameter  of  0.85  pm  to  model  both  cases.  For 
convenience,  we  kept  this  value  constant  for  the  initial  calculations,  rather  than  averaging  over  a 
diameter  distribution.  Thus,  although  we  label  the  graphs  B/r  or  K12,  we  are  primarily 
determining  how  a  substantial  change  in  average  length  affects  D(Q,A )  for  cells  that  are 
otherwise  identical. 

After  the  calculations  for  a  single  width,  we  compare  the  results  with  those 
obtained  over  a  distribution  of  widths  and  lengths  that  are  assumed  uncorrelated  with  each  other. 
We  used  a  width  distribution  similar  in  shape,  but  narrower  than  the  experimental  distribution 
given  in  the  left  hand  column  of  Table  3.  For  both  E.  coli  B/r  and  the  K12  variant,  we  used  the 
same  width  distribution  as  given  in  the  middle  column  of  Table  3.  For  E.  coli  B/r,  we  used  the 
length  distribution  summarized  in  Table  2,  while  for  K12,  we  used  the  length  distribution  given 
in  Table  4. 


19 


Table  2.  Length  Distribution  for  E.  coli 

B/r. 


EcoS  B/r 


40 


length  (jut< 


Figure  4.  Length  distribution  for  E.  coli 
B/r. 


E.  coli  K12 


40 
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Figure  5.  Length  distribution  for  E.  coli 
K12. 


Length 

(pm) 

Percent 

1.75 

16 

2.25 

37 

2.75 

27 

3.25 

16 

3.75 

4 

Total 

100 

Table  3.  Diameter  Distribution. 
E.  coli  B/r. 


Experimental 

Diameters 

(pm) 

Diameters 
Used  in 
Calculation 
when 

Averaging 

Percent 

0.85 

0.7 

14 

0.9 

0.75 

36 

0.95 

0.8 

36 

1.00 

0.85 

14 

Table  4.  Length  Distribution 
for  E.  coli  K12. 


Length 

(pm) 

Percent 

4.25 

17 

5.75 

37 

7.25 

32 

8.75 

11 

10.25 

3 

Total 

too 

The  calculations  were  performed  by  the  same  procedures  as  those  used  for  B. 
cereus  using  only  the  hemispherically  capped  cylinder  model.  For  each  calculation,  we  averaged 
over  the  length  distributions,  and  also  averaged  over  orientation  using  a  9  x  15  mesh  in  #and  </>, 
respectively. 
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We  did  some  calculations  using  the  somewhat  artificial  distribution  of  diameters 
in  the  middle  column  of  Table  3  so  that  we  could  get  good  convergence  to  solutions  when  doing 
a  full  averaging  over  all  lengths  (Table  2)  with  all  diameters  in  the  proportions  indicated  (i.e., 
without  correlation  between  length  and  width).  We  call  this  “£.  coli  B/r  Distribution  (A).”  For 
good  convergence,  the  ratio  of  distance  between  dipole  neighbors  to  wavelength  should  be 
small.6'7  We  compared  the  resulting  graphs  for  D(0,T)  with  the  graphs  for  “£.  coli  B/r 
Distribution  (2/3  A).”  For  the  latter  case,  the  same  diameter  distribution  was  used  with  a  length 
distribution  in  which  each  contribution  to  the  average  had  the  same  weighting  as  that  given  in 
Table  2.  However,  each  length  was  shortened  to  2/3  the  corresponding  value  in  that  table. 

3.3  Spore  Clusters 

A  real  aerosol  of  spores  will  contain  some  clusters  of  spores  as  well  as  single 
spore  particles.  We  modeled  a  small  cluster  of  approximately  5  to  8  spheres  representing  a 
cluster  of  Bacillus  endospores,  which  might  be  found  as  individual  particles  in  an  aerosol.  Each 
individual  spore-sphere  representing  one  spore  in  the  cluster  has  a  diameter  d  of  approximately 
1.0  -  1.5  pm,  and  our  calculations  presented  here  are  those  for  d=  1.2  pm.  The  spore  radius  is 
r  =  dl 2.  The  model  spore  has  a  complex  refractive  index  corresponding  to  the  wavelength 
dependent  values  given  in  Table  1. 

In  an  extremely  simplified  case,  the  spore-sphere  could  be  represented  by  a  single 
dipole,  which  corresponds  to  the  first  term  in  an  accurate  multipole  expansion  for  a  sphere.  To 
examine  the  extent  to  which  such  a  simplified  model  could  represent  the  spores  in  a  cluster,  we 
compared  the  calculated  value  of  the  Mueller  matrix  ratio  S34 ISU  for  that  model  with  that  from 
a  cluster  of  spore  spheres,  each  represented  by  numerous  dipoles. 

Approach  1 :  Simplified  DDA  model  of  small  cluster  of  spores. 

First,  choose  as  parameters  the  spore  sphere  radius  r  and  the  number  n  of  spore 
spheres.  Each  individual  spore  sphere  will  be  represented  by  a  collection  of  electric  dipoles.  We 
build  the  model  using  a  simple  cubic  lattice  to  locate  the  spore  spheres  and  also  use  a  finely- 
spaced  simple  cubic  lattice  grid  to  locate  the  individual  dipoles  making  up  each  spore.  A  super¬ 
sphere  of  radius  R  is  centered  at  one  of  the  lattice  points,  and  R  is  chosen  equal  to  a 
representative  radius  of  the  actual  spore  cluster.  The  spacing  of  the  lattice  for  positioning  the 
spore  sphere  is  then  chosen  so  that  a  closely  packed  set  of  n  spore  spheres  each  of  radius  r  will 
be  within  radius  R  of  the  origin.  The  spacing  between  lattice  sites  is  taken  as 


a  =  \^~  O6) 

thereby  distributing  the  available  space  among  the  spores.  However,  to  avoid  overlap  of  spores, 
the  radius  of  each  spore  is  then  assigned  to  be  half  of  the  lattice  spacing,  thereby  assuring  that  the 
spore  spheres  just  barely  touch. 
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A  3-dimensional  grid  is  defined  corresponding  to  a  finely-spaced  simple  cubic 
lattice.  The  procedure  for  constructing  this  grid  is  the  same  as  that  in  our  earlier  DDA 
calculation.2  3  8  The  desired  number  N  of  grid  points,  typically  several  thousand,  to  fit  into  the 
super-sphere  is  specified  as  input.  Then  the  grid  lattice  spacing  is  taken  as 


agnd  ~ 


3 


N 


(17) 


This  gives  approximately,  but  usually  not  exactly,  the  requested  number  of  grid  points  in  the 
super-sphere.  We  have  shown  elsewhere3  that  using  several  thousand  dipoles  to  represent  a 
single  sphere  gives  results  for  S34  / Sn  very  closely  approximating  those  computed  using  the 
exact  Mie  expansion. 

We  use  the  following  procedure  for  the  calculation:  Each  lattice  site  that  falls 
within  the  super-sphere  has  a  spore  sphere  centered  on  it.  The  radius  of  each  spore  sphere  is  half 
the  lattice  separation.  If  the  spore  spheres  are  {s;-;  /  =1, ...  n},  all  points  <  r  from  each  of  the 

center  points  s-  are  inside  the  /  spore  sphere  with  a  dipole  located  on  each  of  them.  All  other 

points  are  in  air  and  are  assumed  unoccupied  by  dipoles. 

The  polarizability  of  the  DDA  dipoles  is  obtained  from  the  Clausius-Mossotti 
relation,  eq.  14,  taking  the  density  to  be  that  of  the  grid  lattice. 

The  scattering  is  then  calculated  as  a  function  of  angle  using  the  DDA  by  solving 

the  equation 


£(/£)  =  £„(/<?)  +  fjfjaT(i^jo)EUu)  (18) 

7=1  u=\ 

as  in  our  earlier  calculations.  Here  E(ju)  is  the  ^component  of  the  induced  electric  field  at 
dipole  site  j,  E(i< !;)  is  the  £  component  of  the  incident  electric  field  at  the  dipole  site  /,  a  is  the 
polarizability,  and  T{i^\ju)  is  the  component  of  the  dipole-dipole  tensor  coupling  between 

site  i  and  site/.  That  is,  the  field  component  £(/£)  produced  by  dipole  components  at 

each  site  j  for  u  =  x,y,z  is 


09) 

J  u= 1 

Approach  2:  Oversimplified  model  of  small  cluster  of  spores 

In  the  DDA  model,  the  dielectric  particle  is  often  thought  of  as  represented  by  a 
collection  of  small  polarizable  spheres,  with  only  the  dipole  moment  of  each  sphere  to  be  kept 
(making  it  equivalent  for  light  scattering  purposes  to  the  corresponding  collection  of  polarizable 
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point  dipoles).  If  we  took  the  imaginary  small  spheres  in  that  construction  to  be  the  spore 
spheres,  and  really  kept  only  the  dipole  moments  in  any  multipole  expansion  of  the  fields,  this 
would  involve  representing  each  spore  by  a  single  dipole  moment. 

Moreover,  the  packing  fraction  of  the  spheres  would  be  0.5.  This  corresponds 
closely  to  the  packing  fraction  measured  for  spore  clusters  [unpublished  experimental  data.  Dr.  J. 
Bottiger,  U.S.  Army  Edgewood  Chemical  Biological  Center  (ECBC),  Aberdeen  Proving  Ground 
(APG),  MD]. 

3.4  Solution  of  the  DDA  Equations 

We  modified  the  procedure  followed  in  our  earlier  work,1 2  which  involved  storing 
the  complete  set  of  dipole-dipole  tensors  for  coupling  possible  pairs  of  dipoles.  Now,  because  the 
grid  may  have  such  a  large  number  of  sites,  most  of  them  unoccupied  by  dipoles,  the  dipole- 
dipole  matrix  for  interactions  between  all  possible  sites  is  too  large  to  store.  Therefore,  when  the 
program  needs  to  calculate  the  dipole-dipole  matrix  contracted  with  some  3/V-dimensional 

N  3 

vector,  such  as^^aT(i^;  j  u)  E(j  u)  .it  goes  through  a  table  of  which  sites  have  dipoles  on 

7=1  u= 1 

them  and  sums  only  over  those.  It  does  not  store  the  dipole-dipole  interaction  between  each  such 
pair  of  sites,  but  calculates  it  afresh  each  time. 

The  program  offers  an  option  of  solving  the  system  of  equations  by  simple 
iteration  or  by  the  conjugate  gradient  method.  The  simple  iteration  procedure  is  usually 
computationally  less  costly  when  it  converges,  but  the  conjugate  gradient  method  usually 
converges  when  the  simple  iteration  method  fails,  typically  for  larger  particles,  larger  refractive 
index,  and  smaller  wavelength. 

3.5  Calculations  for  Spore  Clusters 

In  Section  4.4,  we  present  a  calculation  that  tests  the  idea  of  using  a  single  dipole 
for  each  sphere  compared  to  using  spore  spheres  in  a  symmetrical  arrangement,  and  we  find  that 
a  useful  rough  approximation  can  be  obtained  in  that  way. 

To  construct  more  realistic  arrangements  of  spores  into  a  cluster  than  placing 
them  on  a  cubic  array,  we  devised  computer  algorithms  to  use  the  method  of  Witten  and  Sander  ' 
to  construct  somewhat  random  arrangements  of  spores  forming  the  cluster  in  the  previous  two 
cases  (dipoles  at  precise  lattice  points,  and  DDA  spheres  at  precise  lattice  points).  We  used  the 
following  algorithm: 

(1)  Place  the  spheres  at  random  locations,  each  time  checking  that  the  spheres 
were  not  assigned  to  overlap.  If  any  added  sphere  overlaps  other  spheres,  reject  the  new  location 
and  try  again. 


(2)  For  each  sphere 

(a)  Displace  all  the  spheres  randomly  by  up  to  8  times  its  radius  vector. 
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(i)  If  any  sphere  then  overlaps  other  spheres,  do  not  accept  the  new  locations,  but 
try  again  to  find  a  new  non-overlapping  configuration.  (The  retry  can  occur  up  to  2000  times.) 

(ii)  For  an  accepted  configuration  in  (i),  if  the  total  displacement 

squared  ^|r,|‘  from  the  origin  increases,  do  not  accept  the  new  location  but  try  again  at  step  (i) 

i 

(up  to  500  times) 

(iii)  For  an  accepted  configuration  in  (i),  if  the  total  displacement  squared 
^|i*  |"  decreases,  accept  these  as  the  new  positions.  Then  go  to  step  (i)  to  try  again  for  an  even 

i 

lower  total  displacement  squared. 

(iv)  If  no  further  reduction  in  the  total  displacement  squared  occurs  within  500 
consecutive  trials,  reduce  8  by  multiplying  by  0.99  and  go  to  (b)  below. 

(b)  Test  if  8  is  large  enough  (>1010)  and  restart  the  search  at  step  (a)  if  it  is.  But 
if  8  has  been  reduced  to  less  than  10"10  print  out  the  locations  of  the  spheres  and  exit. 

This  procedure  thereby  uses  a  random  walk  biased  toward  particle  coalescence  to 
construct  the  cluster  of  a  few  spores  (our  version  of  the  Witten  Sander  model).9’10 

The  program  for  treating  the  spores  as  collections  of  DDA  spheres  was  then 
modified  in  two  ways.  First,  we  no  longer  require  the  spore  spheres  to  be  centered  each  on  a  grid 
point.  Instead,  the  output  of  the  program  giving  the  locations  of  the  centers  of  each  spherical 
spore  and  its  radius  (i.e.,  position./)  is  read  in  as  input.  Secondly,  the  spores  are  prevented  from 
overlapping  simply  by  the  test  within  position./ that  causes  it  to  reject  any  configurations  with 
overlapping  spheres. 

The  sites  on  the  grid  that  fall  within  each  assigned  spore  sphere  are  taken  to  have 
a  dipole  on  each  site.  Since  the  spore  spheres  need  not  have  been  taken  to  be  centered  on  grid 
points,  some  may  contain  more  grid  points  than  others,  and  therefore  may  contain  more  dipoles; 
but,  the  number  will  be  approximately  the  same  if  the  sphere  radii  are  considered  equal. 

We  settled  on  two  different  clusters  of  spore  spheres  as  examples  for  determining 
what  effect  differences  in  the  arrangement  of  the  spores  would  have  on  the  graph  of  the 
depolarization.  Graphs  of  the  calculations  are  presented  for  those  clusters  in  Section  4.4.  We  call 
the  first  of  these  the  “Compact  Cluster.”  Its  shape  can  be  seen  in  Figure  6.  The  coordinates  of 
that  cluster  are  given  in  Table  5. 
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Figure  6.  Compact  Cluster.  Five  spheres  form  a  plane  with  one  sphere  on  top  of  the  center. 
The  lowest  sphere  in  the  foreground  is  slightly  below  the  plane  of  the  five  spheres. 


Table  5.  Coordinates  of  Spore-Sphere  Centers  for  Compact  Cluster 


X 

y 

z 

1 

0. 107968 14E- 
03 

-0.75738501E-01 

-0.25377747E-01 

2 

-1.0337081 

0.61300209 

-0.33699078E-01 

3 

1.0404514 

-0.77011471 

-0.25872477E-01 

4 

0.19538628 

0.66890403 

0.95695171 

5 

-0.42222498 

-0.74854324 

-1.0420674 

6 

0.49930960 

0.76658113 

-0.80320375 

7 

0.54819461 

-0.47870546 

1.0590495 

Figure  7  shows  the  other  cluster  for  which  calculations  are  presented,  which  we  refer  to  as  the 
“Loose  Cluster,”  for  which  the  coordinates  are  shown  in  Table  6. 


Figure  7.  Two  views  of  the  Loose  Cluster. 
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Coordinates  of  the  Loose  Cluster  are  presented  below  in  Table  6. 


Table  6.  Coordinates  for  Loose  Cluster 


X 

y 

z 

1 

0. 

0. 

0. 

2 

0.28940 143E-09 

0.33697962E-03 

1.2 

3 

-0.38798669E-06 

.47030620E-07 

-1.2 

4 

-0.334079 12E-02 

-1.1762189 

0.23769309 

5 

0.62789883E-07 

1.2 

-0.41 98207  IE-05 

6 

1.2 

0.181 16370E-03 

0.21654494E-03 

7 

-0.77928239E-1 1 

-0.87341 592E-07 

-2.4 

3.6  Orientation  Averaging 

When  considering  scattering  from  an  aerosol  for  a  non-spherically  symmetric 
cluster  or  individual  particle,  one  must  do  an  orientation  average.  To  include  averaging  over 
orientations,  some  modifications  were  made  to  avoid  needless  recalculations  within  the  loop 
structure  for  the  numerical  integration  over  orientations.  In  addition,  in  earlier  work,  one  of  the 
three  rotations  over  Euler  angles  could  be  taken  to  be  a  rotation  around  a  rotational  symmetry 
axis  of  the  particle,  and  therefore  could  be  omitted.  This  continues  to  apply  to  the  spheroidal  and 
capped-cylinder  shaped  scatterers.  However,  neither  the  Compact  nor  the  Loose  spore  cluster  we 
generated  has  a  symmetry  axis,  so  that  averages  should  be  taken  over  all  three  orientation  angles. 
The  orientation  procedure  of  our  earlier  work3  was  therefore  followed  with  averaging  over  an 
additional  angle  if/  for  rotation  about  an  axis  fixed  in  the  body  of  the  cluster.  The  number  of 
orientations  averaged  over  was  empirically  determined. 

Our  previous  method  solves  the  DDA  equations  at  each  given  angle  #only  for 
two  perpendicular  polarizations,  and  obtains  the  solutions  for  the  various  (j)  at  that  same  6  as 
linear  combinations  of  the  two  solutions,  thereby  greatly  reducing  the  computer  time  and 
allowing  large  mesh  numbers  for  (j).  This  method  was  revised  to  apply  to  the  spore  cluster 
calculation  with  the  third  Euler  angle  ^included.  It  was  not  possible  in  the  case  of  spore  clusters 
to  take  advantage  of  any  additional  inversion  symmetry.  The  inversion  symmetry  allows  the 
DDA  values  to  be  solved  for  each  #only  over  the  range  0  <  6  <  ;r/2  rather  than  0  <  6  <  n  .  A 
nominal  11x21  mesh,  which  we  used  for  disks,  spheroids,  and  spheres  throughout  this  report 
(and  in  our  earlier  work),  with  the  inversion  symmetry  used  in  this  way,  is  equivalent  to  a  21  x 
21  mesh  when  the  symmetry  advantage  is  not  used. 

Figure  8  shows  the  orientation-averaged  graph  of  S34/  Su  for  over  a  7  x  7  x  5 
mesh  in  the  angles  0,  (/),  and  if/,  respectively,  compared  with  the  results  of  averaging  only  over  6, 
(f>,  and  with  the  results  for  a  single  orientation.  The  calculation  is  for  the  Compact  7-spore  cluster 
with  228  to  251  dipoles  in  each  spore  of  radius  0.6  pm  at  wavelength  1.551  pm  and  with 
refractive  index  n  =  1.503  +  0.004  /.  The  average  over  all  three  angles  is  reasonably  close  to 
average  over  only  two  angles  ^>and  6,  but  both  curves  differ  strongly  from  the  scattering 
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pattern  for  a  single  orientation.  The  averaging  over  #and  ^reduces  the  sharpness  of  the  peaks 
relative  to  those  of  the  single  orientation,  as  might  have  been  expected.  A  stronger  reduction 
results  when  the  additional  averaging  over  the  third  Euler  angle  if/  for  rotations  about  the  z  body- 
fixed  axis  of  the  scatterer  is  included.  We  found  that  the  change  as  a  result  of  averaging  is 
stronger  for  smaller  clusters.  Overall,  we  found  that  scattering  for  the  orientation  average  can 
differ  substantially  from  that  of  individual  spore  clusters. 


. Single  orientation 

—  7*7  orientation  mesh 


Figure  8.  Single  random  orientation  vs.  average  for  the  Compact  7  spore-sphere  cluster. 

. 7x7x5  orientation  mesh 

—  •9x21  x  9  orientation  mesh 


Figure  9.  Comparison  of  various  mesh  sizes  for  the  Compact  7  spore-sphere  cluster. 
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Various  mesh  numbers  were  tried  to  seek  convergence,  as  shown  in  Figure  9  for 
spore  clusters  and  in  Figure  10  for  capped  cylinders. 

The  results  of  several  tests  suggest  that  a  9  *  15  mesh  for  #and  <j> ,  respectively,  is 
adequate  for  scatterers  like  the  capped  cylinder  model,  in  which  the  third  angle  is  unnecessary 
because  of  the  symmetry  axis.  Representative  calculations  are  shown  in  Figure  10.  Note  that  the 
graphs  in  that  figure  overlap. 


Figure  10.  Depolarization  ratio  for  capped  cylinder  of  length  1.5  pm  and  radius  0.35  pm 
averaged  only  over  orientation,  at  wavelength  1.551  pm,  n  =  1.503  +  0.000402  i. 


4.  RESULTS  FOR  INDIVIDUAL  PARTICLES 

In  almost  any  aerosol  encountered,  a  very  large  portion  of  the  particles  will  be 
single  particles.  Methods  are  available  to  build  a  single  aerosol  analyzer  that  can  separately 
measure  the  depolarization  vs.  angle  for  each  single  particle  drawn  through  the  instrument. 

4. 1  Spores 

In  this  section,  we  examine  the  graphs  for  depolarization  obtained  from  single 
particles  that  have  parameters  appropriate  for  Bacillus  spores  modeled  in  two  ways.  The  results 
for  a  capped  cylinder  shape  and  a  prolate  spheroid  shape  are  compared  using  the  same  optical 
and  size  parameters.  Size  averaging  is  accomplished  using  the  parameters  described  in  Section  3. 
Orientation  averaging  is  also  accomplished  as  described  in  Section  3.  We  used  approximately 
2000  individual  dipoles  to  model  the  shape  indicated  for  each  calculation. 
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Figure  1 1.  Mueller  matrix  element  Sn  for  an  ensemble  of  spores  averaged  only 
over  orientations  for  the  single  size  spore  with  length  1 .25  pm,  and  diameter 
or  minor  axis  of  0.72  pm.  The  scattering  wavelength  is  266  nm. 

In  Figures  1 1  and  12,  we  show  the  two  Mueller  matrix  elements  Sn  and  S22  used 
to  calculate  the  depolarization  D(0,A)  for  scattering  wavelength  266  nm.  Figure  13  shows  the 
resulting  D(0,/i). 


Figure  12.  Mueller  matrix  element  Sn  for  an  ensemble  of  spores  averaged  only  over  orientations 
for  the  single  size  spore  with  length  1 .25  pm,  and  diameter  or  minor  axis  of  0.72  pm.  The 
scattering  wavelength  is  266  nm. 
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Figure  13.  Depolarization  for  ensemble  of  spores  of  single  size  with  length  1.25  pm  and 
diameter  0.72  pm  averaged  only  over  orientation,  with  no  size  average.  The 
scattering  wavelength  is  266  nm. 


The  results  for  the  capped  cylinder  and  the  spheroid  models  show  similar 
magnitudes;  but,  the  details  of  the  graph  are  quite  different,  particularly  in  the  forward  direction. 

The  set  of  graphs  shown  in  Figure  14  compares  the  depolarization  results  at 
several  different  scattering  wavelengths  for  the  capped  cylinder  vs.  the  prolate  spheroid  model 
when  averages  are  taken  over  the  semi-empirical  size  distribution  as  well  as  over  orientation. 

. D( 266  nm)  -  capped  cylinder 

—  D(266  nm)  -  spheroid 
—  D(410  nm)  -  capped  cylinder 
—  •  D(410  nm)  -  spheroid 


Figure  14.  Depolarization  ratio  comparison  for  two  spore  models  averaged  over 
size  and  orientation  at  266-,  410-,  and  1551-nm  wavelengths. 
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Comparing  the  results  at  266  nm  with  Figure  13,  we  see  that  the  size  averaging 
produces  some  smoothing,  as  would  be  expected,  and  that  including  larger  sizes  increases  the 
depolarization.  In  this  modeling,  we  used  a  single  width  for  each  length  category.  Since  any 
correlation  between  width  and  length  is  small,  a  more  realistic  model  to  fit  experimental  data 
would  average  over  all  diameters  for  each  length.  This  would  produce  a  much  smoother  graph, 
particularly  at  short  wavelengths  (e.g.,  266  nm).  As  the  scattering  wavelength  increases,  we  see 
that  graphs  for  both  shapes  have  less  structure  as  the  wavelength  increases. 

Figure  15  shows  the  size  and  orientation  averaged  depolarization  graphs  for 
scattering  wavelengths  355  and  750  nm. 


. D(355  nm)  -  capped  cylinder 

—  D(355  nm)  -  spheroid 


Figure  15.  Depolarization  ratio  comparison  for  two  spore  models  averaged 
over  size  and  orientation  at  355-  and  750-nm  wavelengths. 

This  may  be  compared  with  the  results  for  460  nm  shown  in  Figure  16. 
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Figure  1 6.  Depolarization  ratio  comparison  for  two  spore  models  averaged 
over  size  and  orientation  at  460-nm  wavelength. 

Finally,  in  Figure  17,  the  long  wavelength  results  for  1551  and  3389  nm  are 

shown. 

. D(1551  nm)  capped  cylinder 

—  D(1551  nm)  spheroid 
—  •  D( 3389  nm)  capped  cylinder 
—  D( 3389  nm)  spheroid 


Figure  17.  Graphs  for  the  depolarization  at  IR  wavelengths  for  averaging  as  in  Figure  14. 
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We  note  that  the  gap  between  the  graphs  for  the  two  models  appears  largest  at 
wavelength  410  nm  as  shown  in  Figure  14.  In  Figure  17,  we  plot  the  graphs  for  the  two  longest 
wavelengths  studied  (1551  and  3389  nm). 

We  note  also  that  the  depolarization  magnitude  continues  to  decrease,  and  the 
structure  becomes  less  detailed  as  the  size  to  wavelength  ratio  continues  to  decrease.  This  is 
expected  for  scattering  functions  in  general.  From  the  above  examples,  for  wavelengths  near 
4 1 0  nm,  we  see  that  the  results  for  the  two  shape  models  are  quite  distinct  even  though  the  size 
and  index  parameters  are  identical.  If  angle  dependent  data  were  available,  it  appears  that  two 
similar,  but  distinctly  shaped,  aerosols  with  identical  size  and  optical  parameters  could  be 
distinguished  using  depolarization. 

The  direct  backscattering  data,  which  is  important  for  remote  detection 
applications,  is  shown  in  Figure  18.  In  this  case,  the  models  for  the  two  shapes  produce  very 
similar  results,  so  it  is  unlikely  that  two  aerosols  that  are  so  similar  in  their  shape  and  optical 
parameters  could  be  distinguished  on  the  basis  of  backscattering  depolarization.  However,  the 
differences  in  the  depolarization  graphs  are  large  at  certain  wavelengths,  so  an  angular  scattering 
graph  does  reveal  shape  information  for  bioparticles  with  similar  parameters. 

O  Spheroid 
•  Capped  Cylinder 


Figure  18.  Backscattering  (0  =  180°)  depolarization  for  two  particle  shape  models  of  single 
spores,  both  models  averaged  for  size  and  orientation  with  B.  cereus  spore  parameters.  From  the 
left,  the  wavelengths  are  266,  355, 410,  750,  1551,  and  3389  nm.  The  backscatter  results  for  both 
models  overlap  at  410  nm,  and  also  at  3389  nm. 
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4.2 


Vegetative  Bacteria 


In  this  section,  the  graphs  of  D{®,X)  are  presented  for  rod-shaped  bacteria  with 
size  and  optical  parameters  modeled  using  parameters  similar,  but  not  identical,  to  those 
measured  experimentally  for  two  different  strains  of  Escherichia  coli  (Section  3.2  and  below). 

The  calculations  summarized  here  predict  different  behaviors  of  the 
depolarization  ratio  as  a  function  of  angle  for  the  two  strains  of  rod-shaped  bacteria.  We  will 
label  the  graphs  as  K12  or  B/r,  corresponding  to  the  length  distribution  appropriate  to  that  strain 
of  E.  coli.  Details  of  the  graphs  shown  in  Figures  19-22  may  need  some  corrections  because 
calculations  to  assure  the  DDA  completely  converged  were  too  time  consuming  to  complete  on 
the  computers  available. 

Figures  19  and  20  show  graphs  for  Mueller  matrix  elements  Su  and  S22  vs.  angle 
for  the  two  strains,  which  differ  only  in  their  length  distributions.  The  length  distributions  are 
given  in  Tables  2  and  4,  and  the  index  used  in  both  cases  was  n  =  1.373  +  0.00097  /  as 
appropriate  for  both  strains  at  visible  wavelengths.  The  K 1 2  bacteria  are  much  longer  than  the 
B/r  strain.  A  single  diameter  of  0.85  pm  was  used  in  both  cases  to  concentrate  on  how  changes  in 
length  distribution  alone  affect  the  graphs  of  Su,  S22,  and  the  depolarization. 


Figure  19.  Mueller  matrix  elements  Sn  and  S22  vs.  angle  at  wavelength  355  nm 
for  rod-shaped  bacteria  with  Escherichia  coli  B/r-like  parameters,  averaged  over 
length  and  orientation.  The  dashed  line  is  for  S22- 
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Figure  20.  Mueller  matrix  elements  S\  \  and  S22  vs.  angle  at  wavelength  355  nm 
for  rod-shaped  bacteria  with  K  12-like  parameters.  The  averages  are  again  over  both 
length  and  orientation. 

Inserting  the  values  for  Su  and  S22  into  the  equation  for  the  depolarization  leads 
to  the  graphs  for  depolarization  vs.  angle  shown  in  Figures  21  and  22. 


Figure  21.  Depolarization  vs.  angle  at  wavelength  355  nm  for  rod-shaped  bacteria 
with  B/r-like  parameters  (shorter  lengths).  The  averages  are  over  lengths  and 
orientations,  with  no  average  over  diameters. 
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Figure  22.  Depolarization  vs.  angle  at  wavelength  355  nm  for  rod-shaped  bacteria 
with  K1 2-like  parameters  (longer  lengths).  The  averages  are  over  lengths  and 
orientations,  with  no  average  over  diameters. 


In  the  present  case,  we  can  see  that  differences  in  length  distribution  alone 
apparently  give  rise  to  qualitative  differences  in  the  corresponding  graphs  of  depolarization  vs. 
angle.  We  found  that  corresponding  graphs  for  400-nm  wavelength  (not  shown)  are  very  similar 
in  their  gross  features  to  the  graphs  shown  for  355-nm  wavelength  except  for  a  shifting  of 
details. 


The  most  obvious  effect  of  the  increase  in  length  at  wavelengths  in  the  355-  to 
400-nm  range  is  that  the  depolarization  by  the  longer  bacteria  is  generally  larger,  and  the 
envelope  of  that  fluctuating  function  is  fairly  constant  from  30°  to  about  160°  as  shown  by  the 
model  calculations  graphed  above  in  Figure  22.  However,  the  envelope  of  the  function  for  the 
shorter  bacteria  is  generally  rising  from  0°  to  about  160°  as  seen  in  Figure  21. 

Figure  23  shows  side-by-side  comparisons  of  the  depolarization  for  750-  and 
1551-nm  scattering  wavelengths.  Again,  more  structure  is  seen  for  the  longer  strain  at  both 
wavelengths,  while  there  is  a  good  deal  of  similarity  between  the  graphs  for  the  shorter  and 
longer  bacteria  at  the  longer  1551-nm  wavelength. 


36 


Figure  23.  Depolarization  vs.  angle  at  longer  wavelengths  for  rod-shaped  vegetative  bacteria. 
B/r-like  parameters  (shorter  bacteria)  are  on  the  left  and  K  12-like  parameters  (longer  bacteria) 
are  on  the  right.  The  averages  are  over  lengths  and  orientations,  with  no  average  over  diameters. 
The  diameter  was  kept  constant  at  0.85  pm. 

Experimentally,  the  lengths  and  diameters  are  uncorrelated,  so  it  is  desirable  to  do 
an  averaging  procedure  where  every  length  is  averaged  with  every  diameter.  The  two  length 
distributions  used  were  E.  coli  B/r  Distribution  (A)  (average  length  2.52  pm  )  and  E.  coli  B/r 
Distribution  (2/3  A)  (average  length  1.68  pm).  These  distributions  are  defined  in  the  paragraph 
after  Table  4  and  are  used  together  with  an  artificially  small  diameter  distribution  (Table  3, 
middle  column),  which  was  used  to  keep  spacing  between  dipoles  small  to  avoid  convergence 
problems. 


The  resulting  graphs  for  D(©,A)  vs.  0  are  shown  in  Figure  24.  The  depolariza¬ 
tion  is  smaller  for  the  shorter  cells.  This  is  to  be  expected  because  the  limiting  case  as  the 
cylindrical  portion  of  the  model  of  a  cell  gets  shorter  is  a  sphere,  and  the  depolarization  for  a 
sphere  is  identically  zero.  This  is  mathematically  true  for  a  perfect  sphere;  however,  this  should 
be  tested  with  a  dipole  sphere  to  see  how  close  to  zero  the  graph  becomes.  There  appear  to  be 
some  qualitative  differences  in  graphs  as  the  distribution  gets  shorter.  A  more  extensive  study  to 
fully  understand  this  is  warranted. 
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Figure  24.  Graphs  of  depolarization  vs.  angle  with  averaging  over  orientation  along  with 
averaging  over  every  width  together  with  every  length  for  distributions  defined 
in  Section  3.1.  The  parameters  differ  only  in  lengths  used. 

4.3  Particles  that  have  Almost  Zero  Depolarization 

It  is  an  important  observation  that  perfectly  spherical  particles  produce  exactly 
zero  depolarization.  This  fact  arises  because  the  S3  and  .S4  elements  of  eq.  2  are  identically  zero 
for  perfect  spheres,  as  may  be  verified  by  considering  the  symmetry  of  the  sphere.  The  extent  to 
which  natural  background  particles  have  small  enough  departure  from  sphericity  to  make  this 
observation  useful  for  interpreting  scattering  studies  of  biological  aerosols  mixed  with  such 
backgrounds  is  quite  an  interesting  question;  but,  it  has  not,  to  our  knowledge,  been  studied.  This 
would  be  an  interesting  subject  for  a  future  investigation.  A  mist  of  very  small  liquid  particles 
would  have  this  property,  which  could  be  of  practical  interest  for  an  aerosol  consisting  of 
biological  particles  mixed  with  liquid  droplets.  If  we  had  such  a  mist  of  perfect  spheres  mixed 
with  the  biological  particles,  the  correct  formula  for  D(0,T)  would  be 


D(0,A) 


^  ^'l  1  >  bacteria  ^  ^22  >  bacteria 


<  5,,  > 


bacteria 


+ 


r\-p' 

~T 


<su  > 


sphere 


(20) 


where  /?  is  the  fraction  of  biological  particles  in  the  mixture. 

4.4  Clusters  of  Spores 

A  limited  number  of  numerical  simulations  were  done  on  models  of  bacterial 
spore  clusters.  First,  we  present  a  test  of  a  rather  oversimplified  model  of  a  spore  cluster  that  may 
give  adequate  results  in  certain  limited  cases.  The  idea  is  that  spore  spheres  may  be  represented 
by  single  dipoles  for  certain  restricted  calculations.  The  rationale  is  as  follows:  The  scattering 
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from  a  collection  of  dielectric  spheres  can  be  precisely  calculated  by  considering  the  electric 
field  from  each  sphere  expressed  as  a  multipole  expansion  and  summing  to  find  the  total  field. 
The  result  can  be  roughly  approximated  by  keeping  only  the  dipole  moment  contribution  from 
each  sphere.  This  is  mathematically  equivalent  to  calculating  the  light  scattering  from  the 
corresponding  collection  of  polarizable  point  dipoles.  We  consider  here  how  good  an 
approximation  results  when  we  replace  each  sphere  by  a  single  dipole  for  a  collection  of  seven 
spore-spheres. 


To  test  this,  we  calculated  the  Mueller  matrix  ratio  S34/S11  vs.  angle  for  a  cluster 
of  seven  spore-spheres  each  modeled  with  approximately  400  dipoles  and  compared  the  resulting 
graph  with  that  for  seven  individual  dipoles.  The  spore-spheres  were  placed  symmetrically  on  a 
simple  cubic  lattice.  The  center  sphere  was  placed  on  a  particular  lattice  point,  and  then  the  six 
other  spore-spheres  were  placed  symmetrically  in  the  six  nearest  neighbor  points  around  the 
center  sphere  and  tangent  to  it  with  simple  cubic  symmetry  (along  the  x,  y,  and  z  axes). 
Orientation  averaging  was  not  used  in  this  comparison.  One  symmetry  axis  of  the  globule 
(through  three  spheres)  was  lined  up  with  the  laser  light.  Each  of  the  added  spheres  was  tangent 
to  the  center  sphere  but  did  not  overlap  it.  The  scattering  calculated  for  the  cluster  of  seven 
spore-spheres  was  then  compared  with  that  for  seven  single  dipoles  each  placed  in  the  exact 
center  point  chosen  for  each  of  the  spore-spheres.  The  polarizability  a  of  each  of  the  single 
dipoles  was  taken  to  be  the  same  as  that  of  the  spore  sphere  it  replaces.  It  was  calculated  in 
Gaussian  units  using  the  following  equation 

,  3"2-l 

a  =  Am  — - . 

n2  +2 

That  is,  we  used  the  expression  for  the  polarizability  in  the  Rayleigh 
approximation  of  a  sphere  of  radius  a  and  index  n  placed  in  a  uniform  electric  field. 
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Figure  25.  Calculation  using  seven  spore-spheres  each  with  approximately  400  dipoles  compared 
with  calculation  using  only  seven  dipoles.  For  the  latter  calculation,  each  spore-sphere  was 
replaced  with  a  single  dipole  at  its  center.  The  parameters  used  are  wavelength  0.6328  pm,  and 
for  the  spores,  n  =  1.201  and  radius  a  =  0.3164  pm.  The  polarizability  of  each  of  the  seven  single 
dipoles  is  that  of  the  spore-sphere  in  a  uniform  electric  field. 
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The  results  of  the  two  calculations  are  shown  in  Figure  25,  based  on  rather 
arbitrary  parameters  chosen  to  represent  the  spore.  The  comparison  shows  that  the  two  results  are 
quite  distinguishable,  but  are  similar  for  the  two  methods.  The  simpler  seven-dipole  calculation 
is  of  course  very  much  faster.  The  example  presented  here  suggests  that  if  model  calculations  for 
larger  collections  of  spores  are  needed,  qualitative  results  for  the  Mueller  matrix  elements  can 
probably  be  obtained  efficiently  with  this  single  dipole  approach  when  the  diameter  of  the 
scattering  spore  is  the  size  of  the  wavelength  or  smaller.  This  hypothesis  would  need  to  be  tested 
for  various  wavelength-to-size  ratios  and  for  a  more  realistic  (i.e.,  higher)  refractive  index  to  get 
an  idea  of  the  range  of  validity  of  this  approach. 

In  Section  3  (specifically.  Tables  5  and  6,  and  Figures  6  and  7),  coordinates  and  a 
picture  are  shown  for  the  Compact  and  for  the  more  extended  Loose  cluster  of  spore-spheres.  We 
would  expect  the  contribution  to  the  depolarization  from  the  individual  spore-spheres  to  be  quite 
small  (even  though  these  are  not  perfect  spheres).  Therefore,  most  of  the  depolarization  will  be 
from  the  macro  arrangement  of  these  spheres  and  not  from  the  individual  spore-spheres.  About 
230  dipoles  were  used  to  form  each  individual  spore-sphere. 

The  depolarization  vs.  angle  for  the  Loose  and  the  Compact  clusters  averaged 
over  the  three  orientation  angles  ( 0,<f),y /)  is  shown  in  Figure  26. 


Figure  26.  Depolarization  vs.  scattering  angle  for  the  Loose  (solid  line)  and  the  Compact  (dashed 
line)  spore-sphere  clusters  of  seven  spores  for  scattering  wavelength  1.551  pm.  The  refractive 
index  of  the  spores  is  taken  as  1 .503  +  0.0004  /.  The  orientation  average  is  done  over  a 
9x21  x  9  grid. 
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From  Figure  26,  we  see  that  D(0,  A)  has  less  structure  when  the  cluster  becomes 
more  compact.  In  real  experimental  aerosols,  almost  all  clusters  formed  appear  spherical  (oral 
communication  with  Dr.  J.  Bottiger,  ECBC,  APG,  MD). 

Figure  27  shows  how  a  substantial  change  in  refractive  index  affects  D(0,A),  by 
comparing  a  Compact  Cluster  consisting  of  spores  with  the  same  structure  consisting  of 
vegetative  particles  (i.e.,  the  refractive  index  is  reduced  ). 

The  dashed  line  is  the  same  graph  appropriate  for  spore-spheres  as  presented  in 
Figure  26.  The  solid  line  uses  the  index  appropriate  for  vegetative  cells.  The  shape  of  the  graph 
appears  similar  but  is  somewhat  attenuated  for  the  lower  index. 

In  Figure  28,  the  graphs  of  D(0,A)for  both  Loose  and  Compact  clusters  are 
shown  for  the  mid-IR  wavelength  3.389  pm. 


Figure  27.  Depolarization  vs.  scattering  angle  for  wavelength  1.551  pm  for  the  Compact  Cluster. 
The  scattering  particles  corresponding  to  the  dashed  line  have  the  index  1.503  +  0.0004  / 
appropriate  to  spores.  The  particles  corresponding  to  the  solid  line  have  the  index  1.36  +  0.00001 
i  appropriate  to  vegetative  bacterial  cells. 
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Figure  28.  Depolarization  vs.  angle  for  Loose  (solid  line)  and  Compact  (dashed  line)  Clusters  for 
scattering  at  wavelength  3.389  pm,  refractive  index  of  spore  spheres  is  1.51 1  +  0.0136  /. 

As  expected,  the  structure  becomes  smoother  at  longer  wavelengths.  In  Figure  29, 
we  show  the  effect  of  a  small  index  change  for  scattering  from  the  Compact  Cluster  at  the  same 
wavelength. 


Figure  29.  Depolarization  vs.  angle  for  scattering  at  wavelength  3.389  pm  from  the  Compact 
Cluster  of  spores.  The  solid  line  is  for  index  1.479  +  0.04  /.  The  dashed  line  is  for  index  1.51 1  + 
0.0136  i  for  the  spore-spheres. 
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4.5 


Scattering  from  Disks 


Circular  disks  of  a  size  that  would  interfere  with  measurements  of  bioaerosols 
were  briefly  considered  as  a  type  of  background  aerosol  particle.  Figure  30  shows  the 
depolarization  ratio  vs.  angle  for  two  different  disk  sizes  and  a  mixture  of  the  two.  The  graphs 
were  orientation  averaged  over  the  two  angles  (6,  0)  using  a  27  x  51  grid.  The  calculations  used 
3871  dipoles  in  each  case. 


5.  CONCLUSIONS  AND  RECOMMENDATIONS 

An  important  result  from  our  calculations  is  the  observation  that  shape  differences 
as  well  as  length  differences  for  bacteria  and  bacterial  spores  can  be  detected  by  measuring  the 
angular  dependence  of  the  depolarization.  It  appears  that  the  depolarization  vs.  wavelength  graph 
for  direct  backscattering  is  not  very  sensitive  to  shape  changes.  Laboratory  experiments  are 
needed  to  verify  how  well  conclusions  based  on  these  models  replicate  measurements  with  real 
bacteria. 


In  studying  two  different  computational  models  for  spores,  longer  wavelength 
generally  gave  rise  to  lesser  depolarization  and  less  structure  in  the  graphs.  This  is  expected 
because  (unless  there  is  large  absorption),  longer  wavelength  is  equivalent  to  smaller  size.  There 
was  a  large  difference  in  graphs  for  wavelengths  around  400  nm.  Further  study  would  be 
required  to  characterize  the  differences  due  to  changes  in  shape. 


Depolarization  ratio  for  small  disk 
with  n  =  1 .45 


Figure  30.  Depolarization  vs.  angle  for  scattering  from  randomly  oriented  disks  using  a  27  x  51 
grid  in  {6,0)  for  the  sizes  shown.  Scattering  wavelength  is  1.064  pm  with  a  refractive  index  of 
n  =  1 .45  appropriate  to  mineral  dust.  The  thickness  of  the  disk  is  h,  the  radius  is  r. 
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In  the  present  studies  of  rod-shaped  bacteria,  the  shorter  cells  produce  less 
depolarization.  This  is  to  be  expected  because  the  limiting  case  as  the  cylindrical  portion  of  the 
model  of  a  cell  gets  shorter  is  a  sphere,  and  the  exact  value  of  the  depolarization  for  a  sphere  is 
identically  zero.  This  is  mathematically  true  only  for  a  perfect  sphere.  So,  this  should  be  tested 
with  a  dipole  sphere  to  see  how  close  to  zero  the  graph  actually  becomes.  There  appear  to  be 
some  other  qualitative  differences  in  graphs  as  the  distribution  gets  shorter.  A  more  extensive 
study  to  fully  understand  this  is  warranted.  Further  studies  with  substantially  more  dipoles  are 
needed  to  determine  whether  adequate  convergence  with  respect  to  dipole  spacing  was  attained 
in  all  of  the  present  numerical  experiments. 

Methods  were  developed  for  modeling  the  depolarization  from  either  clumps  of 
spores  or  bacteria.  It  appears  this  will  be  useful  when  further  explored.  Thus  far,  however,  this 
modeling  has  only  been  applied  to  very  small  clumps  of  spherical  particles.  Considerably  more 
modeling  would  be  required  before  firm  conclusions  are  possible  regarding  what  is  expected 
from  scattering  experiments  with  real  clumps  of  either  spores  or  bacteria. 
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